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Abstract 

We investigate the long distance asymptotics of various correlation functions for the one- 
dimensional spin-1/2 Fermi gas with attractive interactions using the dressed charge formalism. In 
the spin polarized phase, these correlation functions exhibit spatial oscillations with a power-law 
decay whereby their critical exponents are found through conformal field theory. We show that 
spatial oscillations of the leading terms in the pair correlation function and the spin correlation 
function solely depend on A/jp and 2A/cp, respectively. Here A/cp = ir(n^ — n\) denotes the mis- 
match between the Fermi surfaces of spin-up and spin-down fermions. Such spatial modulations 
are characteristics of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. Our key observation is 
that backscattering among the Fermi points of bound pairs and unpaired fermions results in a 
one-dimensional analog of the FFLO state and displays a microscopic origin of the FFLO nature. 
Furthermore, we show that the pair correlation function in momentum space has a peak at the 
point of mismatch between both Fermi surfaces k = Akp, which has recently been observed in 
numerous numerical studies. 

PACS numbers: 03.75.Ss, 03.75.Hh, 02.30.IK, 05.30.Fk 
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I. INTRODUCTION 



Bardeen-Cooper-Schrieffer (BCS) theory was formulated over 50 years ago as a micro- 
scopic theory for superconductivity. One of the ingredients in BCS theory is pairing between 
electrons with opposite momenta and spins, i.e., matching between the Fermi energies of 
spin-up and spin-down electrons. In the phase where the system is partially polarized, Fermi 
energies of spin-up and spin-down electrons become unequal. This leads to a non-standard 



form of pairing which was predicted independently by Fulde and Ferrell Jl], and Larkin and 
Ovchinnikov |2|. Fulde and Ferrell discovered that under a strong external field, supercon- 
ducting electron pairs have nonzero pairing momentum and spin polarization. At about the 
same time, Larkin and Ovchinnikov suggested that the formation of pairs of electrons with 
different momenta, i.e., k and —k + q where q ^ 0, is energetically favored over pairs of 
electrons with opposite momenta, i.e., k and —k, when the separation between Fermi sur- 
faces is sufficiently large. Consequently, the density of spins and the superconducting order 
parameter become periodic functions of the spatial coordinates. This non-conventional su- 
perconducting state is known in literature as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state. 

More recently, theoretical predictions of the existence of an FFLO state in one-dimensional 
(ID) interacting fermions 4] have emerged by employment of various methods, such as 
Bethe ansatz (BA) [ 



density-matrix renormalization group (DMRG) 7Hll|. quantum 



Monte Carlo (QMC) |12J], mean field theory 13Ml6j] and bosonization 17|]. At finite mag- 
netization, it was found by Feiguin and Heidrich-Meisner tJ that pair correlations for the 
attractive Hubbard model in a parabolic trapping potential has a power-law decay of the 
form n pair oc cos^ffloMVM an d the momentum pair distribution has peaks at the mis- 
match of the Fermi surfaces /cfflo = Tt{n\ — n^). Wave numbers for the oscillations were 
numerically found as 7c(n^ — nj for the pair correlation function and as 2ir(n^ — n^) for 
the density difference (n-f — nj.) [8] . The FFLO pairing wave number was also confirmed by 
the occurrence of a peak in the pair momentum distribution corresponding to the difference 
between the Fermi momenta of individual species 0, [ijj]. From mean field theory, it was 
demonstrated that the FFLO phase exists in the large-scale response 
and even for temperatures up to O.lTp 14]. 

On the other hand, critical behavior of ID many-body systems with linear dispersion in 



the vicinities of their Fermi points can be described by conformal field theory. Some time 



ago, the critical behavior of t 
by Bogoliubov and Korepin 



re Hubbard model with attractive interaction was investigated 



18 



2l|. They showed that ID superconductivity occurs when 



the average distance between electron pairs is larger than the average distance between 
individual electrons of these pairs. This means that the correlation function for the single 
particle Green's function decays exponentially, i.e., (V'nsV'i,*) — >" e _n ^ with £ = Vp/A and 
s =t, I, whereas the singlet pair correlation function decays as a power of distance, i.e., 
(Vn^Vra^ijt^i,^) ~~ n ~ 9 ' Here A is the energy gap, and the critical exponents £ and 9 are 
both greater than zero. This criterion is met when the external magnetic field is small, i.e., 
H < H c . Once the external field exceeds the critical value, i.e., H > H c , Cooper pairs are 
destroyed. Thus both of these correlation functions decay as a power of distance and the 
pairs lose their dominance, i.e., electrons become more or less independent of each other. 

So far, theoretical confirmation of the FFLO state in ID still relies on numerical evidence 
of spatial oscillations in the pair correlations. Despite key features of the T = phase 
diagram 



25j for the attractive Fermi gas were experimentally confirmed using finite 



28] , the unambiguous theoretical 



temperature density profiles of trapped fermionic 6 Li atoms 
confirmation and experimental observation of FFLO pairing is still an open problem. As 
remarked in Ref. that the ID FFLO scenario proposed in Ref. jl^] does not apply to 
ID attractive fermions where quantum phase transition from the fully-paired phase into 



the spin polarized 



also see Refs. 



22 



rase does not belong to commensurate-incommensurate university class, 



26). For ID attractive spin-1/2 fermions with polarization jjj, the 



low-energy physics of the homogeneous system is described by a two-component Tomonaga- 
Luttinger liquid (TLL) of bound pairs and excess unpaired fermions in the charge sector 
and ferromagnetic spin-spin interactions in the spin sector 23]. In this paper, we determine 
the critical behavior of the single particle Green's function, pair correlation function and 
spin correlation function within the context of a TLL. We show that the long distance 
asymptotics of various correlation functions provide a microscopic origin of FFLO pairing 
for ID attractive fermions. 

This paper is organized as follows. We derive finite-size corrections for the ground state 
energy of the system in Section [Til In Section IHIl we derive finite-size corrections for low- 
lying excitations and introduce the dressed charge formalism. Integral equations for each 
component of the dressed charge matrix is solved analytically in the strong coupling limit 



|c| 3> 1. In Section HV] we derive correlation functions for different operators and discuss 
the signature of FFLO pairing. Finally, conclusions and remarks are made in Section |VJ 



II. GROUND STATE AND FINITE-SIZE CORRECTIONS 



We consider Nf fermions with SU(2) spin symmetry in a ID system of length L with 
periodic boundary conditions. The Hamiltonian for the spin- 1/2 Fermi gas 3|, |4| is given by 



N 



H 



f Q2 



(1) 



where c < is the attractive interaction strength. This model is one of the most important 
exactly solvable quantum many-body systems. In recent years, it has attracted considerable 



due to evidence of the FFLO state. 



attention from theory [5|, |6|, |22M25| and experiment [28 
Systems exhibiting novel phase transitions at T = are particularly useful in studying TLL 
physics and the nature of the FFLO state. 

The quasimomenta for unpaired fermions and bound pairs are given by kj and A a ± id 
which satisfy the BA equations 



kjL = 2tiI 3 + 2 tan_1 ( 1 1 ) , 
2A a L = 2nJ a + £ 2 tan" 1 (^^) + E 2 tan_1 (^p) 



(2) 
(3) 



where quantum numbers L and J a are given by 



(mod 1), J a = 



iV„ - JV 6 + 1 



(mod 1). 



(4) 



Here d = c/2, and and Nf, denote the number of unpaired fermions and bound pairs, 
respectively. The energy and momentum for this system reads 



A't, 



j=l a=l j=l a=l 



(5) 



We define monotonic increasing counting functions z^(kj) := Ij/L and z^{K a ) := J a j L 
and re-label the variables k — > k u , A — >■ k b , Ij — >■ 7 u j and J Q — >• I ba so that we can express 



the root densities in a general form as 



N b 



u ~ a=l 

d 1 1 Nu I "" 

Pbih) ■= -rrZb(h) = 7- Y] ai{h - K,j) ~ t 5^ a 2( fc & - fc < 

afcb 77 ' — ' ' — ' 

where a n (fc) is defined by 



i=l /3=1 



n c 



(6) 
(7) 

(8) 



7r [rid) 2 + A; 2 

Here /c aj - (for j = 1,2,..., N a and a = u, b) denote the BA roots for unpaired fermions and 
bound pairs in the ground state. 

Using the Euler-Maclaurin formula for contributions up to 0(1/L 2 ) when L > 1, the 
finite-size corrections to the root densities can be written in the generic form as 



v— \ f p 

P*(k a ) = P { a } {k a ) +J2 K ^(ka ~ kp) P L p (kp)dkp 

£j l J —Op, 

— y 



P=u,b 



24L 2 



P=u,b 



pW 



where 



pu {k u ) 
Pt\h 




K(k) 



P L p(-Qp) 



K uu (k) K ub (k) 
K bu (k) K bb (k) 



(a = u, b) 



(9) 



-a x {k) 
-ai(k) —a,2(k) 



(10) 



Here, the Fermi points are denoted by ±Q a . Notice that K(/c) is a symmetric matrix. 
In order to calculate finite-size corrections for the ground state and low energy excitations, 



we introduce the thermodynamic Bethe ansatz (TBA) [29|,|30|, which provides a powerful and 
elegant way to study the thermodynamics of ID integrable systems. It becomes convenient to 
analyze phase transitions and low-lying excitations in the presence of external fields at zero 
temperature. In the thermodynamic limit, the grand partition function is Z = tr(e~ n ^' r ) = 
q-g/t ^ w h ere the Gibbs free energy is given by G = E — HM Z — fin — TS, and is written in 
terms of the magnetization H, the chemical potential \x and the entropy S 30] . Equilibrium 
states satisfy the condition of minimizing the Gibbs free energy with respect to particle 
and hole densities for the charge and spin degrees of freedom (more details are given in 



Refs. 



22 



30H33|). At zero temperature, the ground state properties are determined by the 



dressed energy equations 

-Qf> 



K a/3 {k a - kp)sp(kp)dk fi , {a = u,b), (11) 

where el°'(fc Q ) are given by 





(12) 



ID many-body systems are critical at T = and exhibit not only global scale invariance 
but local scale invariance too, i.e., conformal invariance. The conformal group is infinite 
dimensional and completely determines the conformal dimensions and correlation functions 
when the excitations are gapless 34]. Conformal invariance predicts that the energy per 
unit length has a universal finite-size scaling form that is characterized by the dimensionless 
number C, which is the central charge of the underlying Virasoro algebra 



35 



36l | . From the 



density distributions (jHJ) and dressed energy equations (ITT]) , the finite-size corrections to the 
ground state energy is given by 

^ = c - % E ««. ( 13 ) 

a=u,b 

where (7=1, and v u and Vb are the velocities of unpaired fermions and bound pairs, 
respectively. They are defined as 



I ds a ( k c 



dp a (k D 



<(±<w , (a=o>i)i (14) 



~-±Q a " Pa(<5«) 2TT Pa (±Q c 



where prime denotes the derivative with respect to k a and p a (k a ) = lim^oo 2nz^(k a ). The 
term represents the ground state energy in the thermodynamic limit, i.e., N,L — >■ oo. In 
the strong coupling limit, exact expressions for the velocities can be found in Refs. [3, \$\ - 



III. LOW-LYING EXCITATIONS AND DRESSED CHARGE EQUATIONS 

Critical phenomena of critical systems are described by finite-size corrections for their 
low-lying excitations. The method we use to study correlation functions of the s pin - 1/2 
Fermi gas with attractive interaction follows closely the method set out in Refs. 4CH43]. 



The conformal dimensions of two-point correlation functions can be calculated from the 



elements of the dressed charge matrix Z. Long distance asymptotics of various correlation 
functions are then examined through the dressed charge formalism at the T = 0. Three 
types of low-lying excitations are considered in the calculations of finite-size corrections. 

Type 1 excitation is characterized by moving a particle close to the right or left Fermi 
points outside the Fermi sea. It is equivalent to changing the quantum numbers I a j close 
to 1^ for unpaired fermions (a = u) and bound pairs (a — b). 1^ characterize the Fermi 
points of each Fermi sea and are given by J + = /™ ax + 1/2 and I~ = /™ in — 1/2. The change 
in total momentum from Type 1 excitations is 

AP=^ (15) 

a=u,b 



and the change in energy is 



2, ^ ^IQ*) 



2tt 
~L 



a=u,b ' 
a=u,b 

Here iV + > (N~ > 0) stems from the change in distribution of quantum numbers close 
to the right (left) Fermi points. This type of excitation is commonly known as particle-hole 
excitation. 

Type 2 excitation arises from the change in total number of unpaired fermions or bound 
pairs. It is characterized by the change in quantum numbers 

N a = I+-I-, (a = «,6), (17) 

i.e., AN a = N™ cited - jVs round . 

On the other hand, Type 3 excitation is caused by moving a particle from the left Fermi 
point to the right Fermi point and vice versa. This type of excitation is also known as 
backscattering. It is characterized by the quantum numbers 

&D a = I "Y a , (a = u,b), (18) 

while leaving AN a unchanged. 

All three types of excitations can be unified in the following form of the finite-size cor- 



rections for the energy and total momentum of the system 



AE = T\ \ i ( AiV )*( Z_1 ) VZ_lAiV + '(AZ^ZV'ZAD + v «( N a + N a) J 

\ a=u,b / 

AP = ^ I 'AN AD + iV u AD u + N b AD b + £ Wa (iV+ - JV") J . 

\ a=u,b / 



(19) 
(20) 



Here we use the notations 

f / AD 




AN = , AD 



u 

AD b 



I v u ^ 











y_ , - 1 2—1 ^ uu ^ u ^ Z ub (Q b ) . ^ _ 

(Q«) z bb (Q b ) 

The dressed charge equations are a set of four coupled integral equations that read 

Z uu {k) = 1 - ai(k - \)Z ub {\)d\ (22) 

= ~ I " a x {k - \)Z uu {\)d\ - [ " a 2 (k - \)Z ub (\)d\, (23) 



ib 



Z bu (k) = - / a x {k - X)Z bb (X)d\, (24) 
J-Qb 

/Qu rQb 
a 1 {k-X)Z bu {X)dX- / a 2 {k- X)Z bb (X)dX. (25) 
J— Qb 

Quantum numbers AD U and AD b (fl8|) are chosen based on the conditions given in Eq. (j4]) 
and also on the conditions that AD U = AN u /2 (mod 1) and AD b = AN b /2 (mod 1). 
Combining both conditions together with the definition given in Eq. (fl8|) yields 

An _ AN U + AN b AN U 

AD U = (mod 1), AD b = — — (mod 1). (26) 

When the external magnetic field H is smaller than the critical field, spin excitations for 
this model are gapped. Once H exceeds this critical field, spin excitations become gapless 
and the system becomes conformally invariant. In this spin polarized phase, spin degrees of 
freedom are suppressed due to the ferromagnetic nature of excess unpaired fermions under 
a magnetic field. Therefore, bound pairs and excess unpaired fermions form two Fermi seas 
which can be described by a two-component TLL at low temperatures. Hence conformal 
invariance results in a universal finite-size scaling form of the energy shown in Eqs. (fl3|) 
and ffl9l . and a universal form of the critical exponents of two-point correlation functions 
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between primary fields (0'(x, t)0(x', t')) which are determined by the finite-size corrections 
of the model. Multi-point correlation functions can be derived by taking the product of 
two-point correlation functions. 

When T = 0, the correlation functions of ID systems decay as the power of distance, but 



when T > they decay exponentially. Following the standard calculations in Ref. 43J, the 
conformal dimensions are given by 



2A± = [z uu AD u + Z bu AD b ± ZbbAN ;~J z UbANb ) + 2N?, (27) 

2A± = (z ub AD u + Z bb AD b ± ^^f^ ) ' + 2JVf, (28) 

where (a = u, b) characterize the descendent fields from the primary fields. General 
two-point correlation functions at T = take the form 

(o(x, t)0(o, o)> = ^(-atw + W) (29) 

(x - iv u t) 2A i{x + iv u t)^ (x - iv b t) 2A t(x + iv b t) 2A * 
The exponential oscillating term in the asymptotic behavior comes from Type 3 excitations, 
i.e., backscattering. Quantum numbers for the low-lying excitations completely determine 
the nature of the asymptotic behavior of these correlations. Here we are only concerned 
with the T = case. 

The four dressed charge equations can be broken up into sets of two pairs. Eqs. (122"]) and 
( |23l) constitute one pair, whilst Eqs. ( 12^1) and ( 1251) make up the other. Since we are interested 
in the strong coupling limit |c| ^> 1, both sets of equations can be solved iteratively up to 
accuracy l/\c\. Let us consider the first set. Substituting Eq. (122]) into Eq. (123]) and iterating 
the terms give 

/Qu rQb rQu 

d\ ai (k-\)+ / d\ d\'a 2 (k - A)a 2 (A - A') 
■Qu. J —Qh J — Qu. 



-Qu 

d\ I d\' I d\"a 1 (k - A)a x (A - A')ai(A' - A") + . . . (30) 



-Qu 

The functions a n (k) have leading order l/|c|, hence we can ignore all terms that have two 
or more multiples of a n (k). This procedure yields 

Z ub (Q b ) w - / d\ ai (Q b - A) « — 

J-Qu *\c\ 
Substituting Eq. into Eq. ([22]), we obtain 

/Qb rQu 
dX dX'a 1 {Q u -X)a 1 {X-X') + ... (31) 
■Q b J-Qu 

w 1 (32) 
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FIG. 1. These figures show a plot of the dressed charges Z UU (Q U ), Z u b(Qb), Zb u {Qu) and Zbb{Qb) 
versus polarization for different values of [ —y | . 

Next, we consider the second set of equations. Repeating the same arguments as before, 
Eq. (|25|) at the Fermi point Qb becomes 



ZbbiQb) — i 



d\a 2 {Q b - A) + / dX J dX'a^Qb- X)a 1 (X - A') 

Qb 



dX / dX'a 2 (Q b - X)a 2 (X- X') + ... 

■Q b J-Qb 
1 2Q fe 



7T C 



(33) 



Eq. (l24l) at the Fermi point Q u then reads 



ZbuiQv) 



dXai{Q u - A) + 



dX 



dX'ai(Q u - A)a 2 (A - A') + . . . 



4Q 6 



7T C 



(34) 
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From Ref. 22], the Fermi points in the strongly attractive limit are given by 



(35) 



Q u nirn f P ( 1+ 2 ^| 7 | jP J ' 

where rif — Nf/L is the density of fermions per unit length, 7 = c/n/ is the dimensionless 
interaction parameter and P = (N^ — Ni)/Nf = N u /Nf is the polarization. Inserting these 
relations into the expressions for dressed charges, we obtain 

4P 



171 

fl-F) (1-P) 

z bu (Qu) « -^r^, « 1 - ^r^- ( 3? ) 

|7I 2 l7l 
In FIG. [TJ the dressed charges are numerically calculated and plotted against polarization 
for different values of interaction strength |7|. 

In the strong coupling limit, the external magnetic field H is related to the polarization 



as 



rr ri l 1 1r,2 4(1 - P) AP \ 7T 71 (1 — P) ( AP\ , . 

H « + 2n 2 n 2 P 2 1 + v , , ; - — i '- 1 + — . 38 

2 v i7i m) 8 v 

With this relation, we can evaluate the dressed charges for different values of H. From the 
expressions for the dressed charges in Eq. (|37|) . the conformal dimensions in terms of 
polarization are given by 

2A J „ ( AA , ± ^ 2 _ j£ ( AB „ ± (AD, T ^) + 2^, (39, 



IV. CORRELATION FUNCTIONS AT ZERO TEMPERATURE 

Here we consider 4 types of correlation functions, namely the single particle Green's 
function G^(x,t), charge density correlation function G nn (x,t), spin correlation function 
G z (x,t), and pair correlation function G p (x,t). Each correlation function is derived based 
on the choice of AN U and AN^. 



11 



The one particle Green's function, which is also called the Fermi-field (FF) correlation 
function in some literature, decays exponentially when the external magnetic field is not 
strong enough to overcome the gap associated with the breaking of bound states 18H21]. 
Once in the gapless phase, i.e., when if cl < H < H c2 where H c i and H c2 are the critical 
fields mentioned in Ref. 122J], every correlation function at zero temperature decays spatially 



as some form of power law 



34 



36 



38 



39j. G t (x,t) is characterized by (AN U , AN b ) = (1,0) 



which in turn allows quantum numbers AD U G Z + 1/2 and AD b G Z + 1/2. The leading 
terms are then given by 



G t (x,t) = (^|(x,t)^ t (0,0)) 



A t;1 cos(7r(n t - 2n l )x) Af, 2 cos (vm+x) 

\x + \v u t\ ei \x + iv b t\ 62 \x + iv u t\ 0s \x + iv b t\ 9 ^ 



where the critical exponents are given by 



AP „ 1 (1 - P) AP 

|7| ^ ^171 |7| 

AP 1 (1 - P) AP 

7 2 2 7 7 



The first term in G^(x,t) comes from (AD u ,ADb) = (1/2,-1/2) and the second term 
comes from (AD U , AD b ) = (1/2, 1/2). The constants Af t i and Af t 2 cannot be derived from 
the finite-size corrections for low-lying excitations. Here we only aim to evaluate the long 
distance asymptotics of these correlation functions. Instead of using N u and N b in the 
oscillation term, we choose to use = / L and = N^/L to elucidate the imbalance in 
the densities of spin-up and spin-down fermions. Both sets of variables are related by the 
relations N u = N t - N± and N s = N±. 

Next we consider the charge density correlation function G nn (x,t) together with the spin 
correlation function G z (x,t). Both of these correlation functions are characterized by the 
set of quantum numbers (AN U , AN b ) = (0, 0) which allows quantum numbers AD U G Z and 
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AD b G Z. The leading terms are given by 



G nn (x,t) = (n(x,t)n(0,0)) 



n 2 A nn! iCOs(27r(n t -wJrr) A w » i2 cos (27177,40;) 

\x + i^t] 611 \x + iv b t\ 92 

A raTOi3 cos(27r(n t - 2n x )x) 
\x + iv u t\ 6 3\x + iv b t\ e * ' 



G 2 (a;,t) = (S z (x,t)S z (0,0)) 

, _ N9 A 2 1 COS I-,, , , ,. ,.. , 

~ ("i ) + 1 . ' z — — + 



. { co s (27r(w t - 7i4,)x) cos (27rn ; a:) 

x + iw u t| ei |x + ii^l 6 * 2 



A 2;3 cos (27r(n t - 2n^)x) 
\x + iv u t\°*\x + iv b t\ * ' 1 ' 

where the operators n(x,t) and S z (x,t) are given in terms of the fields as 

n(x, t) = 4>\{x, t)iJ t (x, t) + 4>\(x, t)^(x, t), (45) 
S z (x,t) = X - (^\{x,t)^(x,t)-^l(x,t)^(x,t)^ . (46) 

The critical exponents for asymptotic expressions of G nn (x,t) and G z (x, t) are 

2(1 -P) 



0! « 2, 6*2 « 2 



7 



# 3 ~ 2 + -p-j-, # 4 « 2 - v ; + —. 47 

ItI Iti Iti 

The constant terms for G nn (x,t) and G z (x,t) come from the choice of quantum numbers 
(AD U , AD b ) = (0, 0). The second, third and fourth terms arise from the choices (1, 0), (0, 1) 
and (—1,1), respectively 

Finally we consider the pair correlation function G p (x,t). This correlation function is 
characterized by the set of quantum numbers (AN U , AN b ) = (0, 1) which allows quantum 
numbers AD U e Z + 1/2 and AD b e Z. The leading terms are 

G p (x,t) = (i(,\(x,t)4>l(x, t)^(0, 0)^(0,0)) 

_ A p ,i cos (?r(n t - nj>) A p , 2 cos (?r(n t - 3nJ:r) g 
~ \x + iv u t\ 9l \x + iv b t\ 9 i \x + ivj\^\x + iv b t\ e * ' 1 J 

where the critical exponents are given by 

1 2' 2 2 2| 7 | ' 

1 8P . 5 5(1 -P) 8P 

2 |7| 2 2|7| |7| 
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The first term in G p (x, t) arises from the choice of quantum numbers (AD U , AD^) = (1/2, 0), 
whilst the second term arises from the choice (AD U , AD^) = (1/2,-1). 

The leading order for the long distance asymptotics of the pair correlation function 
G p (x,t) oscillates with wave number Akp, where Akp = 7r(n^ — n±). Meanwhile, the leading 
order for the spin correlation function G z (x,t), which can also be thought of as the corre- 
lation of the density difference between spin-up and spin-down fermions, oscillates twice as 
fast with wave number 2Akp. The oscillations in G p (x,t) and G z (x,t) are caused by an 
imbalance in the densities of spin-up and spin-down fermions, i.e., — nj., which gives rise 
to a mismatch in Fermi surfaces between both species of fermions. These spatial oscillations 
share a similar signature as the Larkin-Ovchinikov (LO) pairing phase Q]. Our findings of 
the wave numbers agree with those discovered through DMRG ffl, QMC Q and mean 
field theory 14j. Though from conformal field theory, we see clearly that the spatial oscil- 
lation terms in the pair and spin correlations are a consequence of Type 3 excitations, i.e., 
backscattering for bound pairs and unpaired fermions. A comparison between our results 



12] suggest that the coefficient A 



P ,i 



and the results from numerical methods in Refs. [7h 
is very much larger than the coefficient A Pt2 because the frequency of the oscillations in 
numerical studies of G p (x,t) is almost identical to ir(n^ — n±). This observation also applies 
to G z (x,t), where A z .2 and A Zt $ are much smaller when compared with A z ^. 

The correlation functions in momentum space can be derived by taking the Fourier trans- 
form of their counterparts in position space. From Refs. 
equal-time correlation functions of the form 



42 



43), the Fourier transform of 



g(x,t = + ) 



exp(ikox) 



x-iO) 2A+ (x + iO) 



2A- 



(50) 



where A ± 



+ A^ is given by 



g(k « fc„) ~ [sign(A; - k )] 2s \k - k, 



(51) 



The conformal spin of the operator is s = A + — A and the exponent v is expressed in terms 
of the conformal dimensions as v = 2(A + + A~) — 1. 

Hence the equal time correlation functions near the singularities fco for the one particle 
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FIG. 2. (Color online) This figure shows a plot of the pair correlation function in momentum 
space G p {k) against k for different values of polarization P when I7I = 10 and total linear density 
nt = 1. The location of the peaks are at k = 0, 0.2-7T, 0.5-7T and 0.87T when P = 0, 0.2, 0.5 and 0.8, 
respectively. 



Green's function, charge density, spin and bound pairs are 

G t (k) ~ [sign(A; - vr(n t - 2n x ))} 2 ^\k - vr(n t - 2 Ui )\ v \ 
G nn (A;) ~ [sign(A; - 2vr(n t - n ; ))] 2s ""|A; - 27r(n t - n^" 
G z (k) ~ [sign(fc - 27r(n t - nj)] 2s *|fc - 27r(n t - njl"', 
G p (fc) ~ [sign(fc-7r(n t -nJ)] 2 ^|A;-7r(n t -nJ|^, 



where the exponents are given by 

4P (1 - P) 



2s t « 1 + 



7 



2Sr 



2s 2 « 0, 



2s p « 0, 



7/ Ri 



171 

(1-P) 



77^ 



1 8P (1 - P) 

2 + ^" 2| 7 | 



(52) 
(53) 
(54) 
(55) 

(56) 
(57) 
(58) 



We would like to stress that the momentum space correlation functions derived in Eqs. (|52|) - 
( 155]) are only accurate when the momenta k are within the proximity of the wave numbers 



15 



ko, i.e., when k ~ k . FIG.|2]plots G p (k) against k as polarization P varies between to 0.8 



121 ] . We stress again 



This figure is in qualitative agreement with the ones given in Refs. {7,9, 
that our plot is accurate only within the vicinity of the singularity, i.e., when k approaches 
7r(n^ — Tij,). We plotted G p (k) for the entire domain k G (0, 7r) so that readers can visualize 
the curves more easily. 



V. CONCLUSION 



In conclusion, we investigated various zero-temperature correlation functions for the spin- 
1/2 Fermi gas with attractive interaction. We derived the finite-size corrections for ground 
state and low-lying excitations of the model. Using conformal field theory, critical exponents 
of the correlation functions were given in terms of polarization and interaction strength. 
We found that the leading terms of the pair correlation function and the spin correlation 
function oscillate with frequencies ix(n^ — n^) and 27r(n-|- — n±), respectively. We also found 
that backscattering between the Fermi points of bound pairs and unpaired fermions results 
in a ID analog of the FFLO state and displays a microscopic origin of the FFLO nature. 
Furthermore, we showed that there is a peak in the pair correlation function in momentum 
space at k = Tr(n^ — n±) which confirms the oscillation frequency. 

In the spin polarized phase, these correlation functions exhibit spatial oscillations with a 
power-law decay. This critical behaviour can be viewed as an analogy to long range order in 
ID, i.e., the power law decay of the pair correlation function which is regarded as evidence 
of a superconducting/ superfluid state. We also like to mention that from the dressed charge 
formalism, the asymptotic behavior of the correlation functions derived in this paper can 
be numerically obtained with high accuracy for arbitrary interaction strength. Addition- 
ally, by considering weakly perturbed inter-tube interactions or inter-lattice interactions (ID 
fermionic Hubbard model), quasi-lD correlations in the spin polarized phase can be calcu- 
lated from perturbation theory |l9] . This provides a promising opportunity to estimate the 
critical temperature for high-Tc superconductors/superfluids by studying ID to 3D trapped 
cold atoms. 
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